Ab-initio electron transport calculations of carbon based string structures 
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First-principles calculations show that monatomic strings of carbon have high cohesive energy and 
axial strength, and exhibit stability even at high temperatures. Due to their flexibility and reactivity, 
carbon chains are suitable for structural and chemical functionalizations; they form also stable ring, 
helix, grid and network structures. Analysis of electronic conductance of various infinite, finite 
and doped string structures reveal fundamental and technologically interesting features. Changes 
in doping and geometry give rise to dramatic variations in conductance. In even-numbered linear 
chains strain induces substantial decrease of conductance. The double covalent bonding of carbon 
atoms underlies their unusual chemical, mechanical and transport properties. 
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PACS numbers: 73.63.-b, 73.22.-f, 73.40.Jn 

While the research on molecular electronics is in 
progress since 1960s, the problem of interconnecting 
those molecular devices has remained as a real challenge. 
Recently, linear gold chains suspended between two gold 
electrodes have been produced |lj . Another ultimate one- 
dimensional nanowire, namely monatomic carbon linear 
chain (C-LC) has been observed at the center of multi- 
wall carbon nanotubes 0- Some carbon chain struc- 
tures have been subjects of earlier theoretical studies 



This Letter predicts new stable structures of carbon- 
based strings and their unusual electronic transport prop- 
erties. Carbon strings having impressive mechanical and 
electronic properties can form helix structures through 
bending, and also networks or 2D grids through T and 
cross-bar junctions. Interesting and potentially useful 
conductance variations can be achieved in these struc- 
tures. In spite of the fact that the parent diamond crystal 
is a good insulator, an ideal C-LC is a better conductor 
than gold chain. Double-bond formation between adja- 
cent C atoms underlies all these unusual properties of the 
carbon chain. 

We carried out total energy and electronic structure 
calculations using first-principles pseudopotential plane 
wave method Q within density functional theory (DFT) 
and supercell geometries. Generalized gradient approx- 
imation (GGA) and local density approximation (LDA) 
have been used, and all the atomic positions and lat- 
tice parameters of string structures have been optimized. 
The analysis of quantum ballistic conductance has been 
performed using transiesta-C, a recently developed ab- 
initio transport software based on localized basis sets, 
DFT and non-equilibrium Green's function (NEGF) for- 
malism [To| . 

Among a large number of nanowire structures of car- 
bon including planar zigzag, dumbbell, and triangular 
structures, C-LC (known as cumulene) is found to be 
the only stable structure 0. Upon relaxation, all other 
structures are transformed to C-LC with a bond length 



of c = 1.27 A. The C-LC is energetically very favor- 
able as reflected by the cohesive energy of as large as 8.6 
eV/atom (that is close to the calculated GGA cohesive 
energy in diamond structure, 9.5 eV/atom). 

Structural and electronic properties of C-LC are es- 
sential for studying its transport characteristics (see Fig. 
0. The cr-band below the Fermi level, Ep, is composed 
from 2s and 2p z atomic orbitals. The doubly degenerate 
7r-band formed from the bonding combinations of 2p x and 
2p y atomic orbitals U, crosses Ep at k = n/2c. There- 
fore the C-LC structure is vulnerable to Peierls distortion 
which may induce a metal-insulator transition. The exis- 
tence and degree of such a distortion has been a subject 
of debate ||. We calculated that the maximum energy 
gain (per atom) through the dimerization of the chain is 
just 2.7 (3.2) meV using LDA (GGA) at bond alterna- 
tion parameter 6 = 0.017A. However, treating the GGA- 
calculated energy-change profile as an effective potential 
for the dynamics of a C atom in the chain, the zero-point 
energy of the atom in the double-well potential turns out 
to be much higher (51 meV) than the Peierls distorted 
ground state energy (see Fig. ^b)). As a result, the 
small distortion effect is insignificant and cannot be ob- 
served even at T=0 K as calculated with GGA and LDA 
energies [l2j . 

In addition to the cr-bond derived from the states of the 
CT-band, half filled 7r-band states form a 7r-bond as shown 
in Fig. ^c). This double-bond structure is responsible 
for the high stability and metallicity of ideal C-LC. Cal- 
culated longitudinal vibrational modes have higher fre- 
quencies than the doubly degenerate transverse modes as 
shown in Fig. C^d). This indicates that the chain is stiff 
along its axis but flexible in the transverse directions. 
Accordingly, C-LC is suitable for deformations to form 
other string structures. The elastic response of the C- 
LC to different deformations is summarized in Fig. ^e). 
As a measure of the elastic stiffness of the C-LC the sec- 
ond derivative of the strain energy per atom with respect 
to the axial strain, d 2 E/de 2 , is calculated to be 119 eV. 
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FIG. 1: (a) Energy band structure of C-LC. Zone-folded 
bands upon doubling the cell size from c to 2c are shown 
by dashed lines, (b) Variation of LDA and GGA total ener- 
gies per unit cell of C-LC in the presence of Peierls distortion. 
The inset shows the zero-point oscillation energy of a C atom 
in the effective potential due to Peierls distortion, (c) Con- 
tour plots of total, cr-band, and 7r-band charge-densities in the 
plane passing through C-LC. (d) Longitudinal (L) and doubly 
degenerate transversal (T) branches of phonons. (e) Varia- 
tion of total energy with the lateral (solid line), transversal 
(dashed line) displacements of a single C atom in a supercell 
of four atoms. Dot-dashed line (U) shows the variation of the 
total energy with uniform extension or compression of c from 
equilibrium. 



This value corresponds to a very high axial strength for 
the C-LC when compared to the corresponding values for 
carbon nanotubes which are in the range 52-60 eV 0] . 

Even more interesting are structural and chemical 
functionalizations of carbon strings as illustrated in Fig. 
12 Carbon rings, C-Rjv, containing N atoms have a num- 
ber of interesting features, and is closely related with the 
helix structure [?J. The helix structure by itself, is de- 
noted as C-HX(iV,p), and can be generated from C-LC 
with different radii and pitch lengths p (in A). Here N 
is the number of C atoms in one pitch of the helix. The 
binding energies of helices are slightly smaller than that 
of the C-LC due to strain energy implemented by the cur- 
vature, such as for p = 10 A fixed, is 8.3 and 8.5 eV for 
N = 10 and 16, respectively. Finite helix structure are to 
be stabilized by fixing from both ends against a weak ten- 
dency of transforming to a LC. All C-HX structures we 
studied here have the bond alternation no matter what 
the value of N is, and are stable even at T — 1200 K. We 
believe that C helix structure may be relevant to study 
chiral currents and to fabricate nanosolenoids. Further to 
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FIG. 2: Carbon-chain-based nanostructures: Carbon linear 
chain (C-LC); bundle of several linear chains (C-B) with a 
repulsive inter-chain interaction; ring structure with N atoms 
(C-Rjv); helix of N atoms with pitch length p (C-HX(A,p)); 
T-junction of two chains (C-T); cross-bar structure of two 
chains (C-CB); planar grid structure (C-G); and C-LC func- 
tionalized with H-atoms ((CH^)™). 



circular ring and helix structures, the C-LC can branch 
off to form T and cross junctions. These junctions can 
form, since an additional C atom can easily be attached 
to any carbon atom of C-LC with a binding energy of 8 
eV. For example, a T-junction is produced if a perpendic- 
ular chain develops from that adatom. At the junction 
the double bond structure changes into planar sp 2 -like 
bonding. In the case of cross-junction, two carbon atoms 
are attached to a single atom of C-LC. At the junction a 
3D sp 3 -like bonding develops as shown in Fig. [21 These 
chains can further branch off to produce a grid of (N x M) 
atomic rectangular cells 0], or a network structure. In 
principle, these grids or networks can be functionalized 
to generate 2D and 3D artificial crystals. Since the in- 
terconnects between the nodes are conducting, one can 
expect that carbon chain networks or grids can be useful 
for the integration in the molecular electronics. 

Carbon chains are reactive owing to their double bond 
structure. Oxygen atom can be chemisorbed to a sin- 
gle carbon atom (specified as top site) or form a bridge 
bond between two carbon atoms. Similarly, a single hy- 
drogen atom can easily be adsorbed to a chain atom with 
Eb = 3.6 eV; it is larger than the energy associated with 
the chemisorption of H on single-wall carbon nanotubes 
[l5| |. The (CH 2 )„ structure in Fig. |2]is a wide band-gap 
semiconductor . 

Having revealed the atomic and electronic structures 
we now turn to the prime objective of our study and 
examine how the conductance of these infinite and fi- 
nite C-based string structures vary with geometry, elas- 
tic deformation and doping. As usual, we have used an 
electrode-device-electrode geometry for conductance cal- 
culations. The conductance of the device has been calcu- 
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FIG. 3: Calculated conductance of various carbon chain struc- 
tures in units of the conductance quantum (Go = 2e 2 /h). (a) 
Infinite structures of undistorted C-LC and two helices C- 
HX(7V, p) with different pitch lengths, (b) Finite-sized helices 
in between two C-LC electrodes. Each panel is divided into 
two; the top parts are for a single period of the helix, and the 
bottom parts are for two periods. The Fermi levels are set to 
zero in all systems. 

lated as G(E) = (2e 2 /h)Ti(T 1 G r r 2 G a ) (G^: retarded 
and advanced Green's functions, Ti^'- coupling functions 
to the electrodes) after iteratively solving NEGF and 
DFT equations of the system . In order to match the 
device potential and the surface potential of the semi- 
infinite electrodes, the device regions are defined to con- 
tain some portion of the electrodes. Here, we choose to 
use C-LC as the metallic electrodes. The atoms in the 
device region are subjected to geometry optimizations in 
supercells, keeping the electrode-atom positions fixed. 

Earlier, conductance properties of short C chains have 
been reported. Using jellium electrodes it has been shown 
that the conductance of C atomic chains varies in an os- 
cillatory manner with the number of atoms in the chain 
[lif . Larade et al 01 predicted negative differential re- 
sistance for C chains coupled with Al electrodes. They 
also point out quantitative differences in odd- and even- 
numbered C chains. We investigate effects of strain, 
impurities and adatoms on conductance of C strings. Re- 
sults of ab-initio electron transport calculations of C-LC, 
infinite and finite segments of C-HX with varying N and 
p are summarized in Fig. [3] The equilibrium conduc- 
tance of ideal, undistorted C-LC is 2Go due to the dou- 
bly degenerate 7r-band crossing the Fermi level. In case 
of bond alternation due to Pcicrls distortion there opens 
a gap around Ep. However, according to the present 
GGA results the weak bond alternation should not be 
observable due to zero-point oscillations of the atoms. 
We also note that any defect (such as displacement of 
atoms from their equilibrium positions, etc.) in C-LC 
should lead to the localization of current transporting 
state, which is characterized by the localization length 



FIG. 4: Conductance variations of finite C-LCs modified by; 
(a) strain under tension, (b) substitutional impurity atoms of 
N and P, (c) chemical doping with single C or H atom; single 
O atom at top site (dashed curve is obtained by deleting the 
adsorbed O but by maintaining the local distortion in C-LC 
caused by O); single O atom adsorbed at the bridge site. The 
Fermi levels are set to zero in all systems. 

£. For a strictly ID wire like C-LC, £ is proportional 
to the elastic mean free path I (t itself depends on the 
type of defect and temperature) . If the length L of the 
C-LC is larger than £, the resistance of the chain be- 
comes R > + ^e L/€ [Hj]. On the other hand the 
opening of a conductance gap around Ep is related to 
the dimerization of carbon atoms, and that effect is seen 
to be reflected in the conductance plots of infinite carbon 
helix structures. The bending of the carbon chain into a 
helix also introduces some narrow conductance drops of 
integer number, which changes with the pitch length and 
radius of the helix. In case of finite-length solenoids, we 
considered single and double periods of helices placed in 
between two semi-infinite C-LC electrodes. The dimer- 
ization of atoms persists for finite-length helices. In Fig. 
|3fb) the set of curves presented for the conductance of 
finite-sized helices all display oscillatory behavior; the fre- 
quency of which depends on the radius and pitch length 
of the helix. Each peak corresponds to a molecular or- 
bital level of the device (helix) aligned with Ep of the 
electrodes (C-LC). As the axial size of the helix is dou- 
bled the number of the peaks in the conductance pattern 
is also doubled, pointing to the limiting case of infinite 
helix structure. 

How the ballistic electron transport depends on strain, 
substitutional impurity and external doping of C-LC is 
another interesting aspect of our study. When an ax- 
ial tension F is applied to an A-atom C-LC the chain 
elongates (the total length scales by 1 + e, where e is the 
strain) resulting with non-uniform bond lengths. For ex- 
ample, for F — 7.78 nN a 6-atom chain attains a strain 
e = 0.15 and a sequence of equilibrium bond lengths; 
1.56, 1.31, 1.59, 1.30, 1.59, 1.31, 1.56 A when relaxed 
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in between two C-LC electrodes. Resulting conductance 
spectrum of the strained chains displays oscillatory be- 
havior with E. Fig. 01a) presents G{E) of a C-LC device 
including N = 5 or N — 6 atoms subjected to differ- 
ent strains. The character of oscillations is dramatically 
different depending on N being odd or even (We have 
confirmed that chains with N = 4 and 7 comply with 
this trend). If N is odd, the chain always has a peak 
conductance value of 2Go at the Fermi level, whereas 
for even N, G(Ep) corresponds to a dip of oscillations 
which drops considerably with increasing e. This salient 
feature is related to the alignment of the molecular en- 
ergy levels of the device region with the electrochemical 
potential of the electrodes. It turns out that Ep coincides 
with the HOMO (highest occupied molecular orbital) for 
odd- numbered chains, and with the center of HOMO 
- LUMO (lowest unoccupied molecular orbital) gap for 
even-numbered ones. When geometrical optimization of 
the chain is omitted and the strain is assumed to be uni- 
form the conductance oscillations are still present but 
their amplitudes are smaller as shown in the insets of Fig. 
0Ja). The observation that even- numbered strained C 
chains have lower conductance values is compatible with 
the findings of Lang and Avouris for unstrained chains 
PH. They report G ^ 1.8G for C chains with N = 5 and 
7 in between jellium electrodes. In our model, use of C- 
LC electrodes naturally corresponds to a "better contact" 
which leads to G — 2Gq. With our choice of electrodes, 
unlike the results in Ref. flil the difference of N being 
odd or even in conductance of unstrained C chains does 
not manifest itself, because for all N the configuration 
is the same. Therefore, the odd - even N disparity we 
find in conductance of C chains is a pure strain-induced 
effect. 

G(E) is also affected by substitutional impurities as 
shown in Fig. Efb). We considered 7-atom C-LCs con- 
taining a substitutional impurity atom of N or P, at the 
center. While the decrease in conductance is minute for 
N-impurity, P causes a substantial decrease. Conduc- 
tance variations G(E) for C-LC externally doped by a 
single (C, H, and O) atom are illustrated in Fig. EIc). 
The adatom deforms the C-LC by changing the type and 
the geometry of the local bonding at the adsorption site. 
For all cases, we see that G(E) dips to zero at Ep. The 
dip is dominantly due to the deformation of the chain 
since it persists even if the adatom is removed but the 
distortion of the chain is kept the same. The G(E) curve 
reflects the local electronic structure at the adatom site, 
and is a fingerprint for the type of the adatom, as well as 
the adsorption geometry. 

In summary, we showed that carbon atoms can form 
stable string structures with impressive physical prop- 
erties. C-LC is flexible but have very high stiffness 
along the axial direction. Accordingly, carbon strings 
are suited for forming various stable nanostructures; they 
can be easily doped and functionalized by chemisorption 



of adatoms and molecules. Calculated conductance is 
found to be sensitive to deformation geometry and dop- 
ing of the chain structure. Strained C chains show an 
odd-even N disparity. As an interconnect carbon atomic 
string structures can be a potential alternative to gold 
chain as well as to metallic carbon nanotubes. We believe 
that with their novel electronic, mechanical, and trans- 
port properties, carbon strings deserve further theoretical 
and experimental research to be a potential material for 
nanoscience and nanotcchnology. 
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